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Abstract: One-loop 0-exact quantum corrections to the neutrino propagator are com- 
puted in noncommutative U*(l) gauge-theory based on Seiberg-Witten maps. Our closed 
form results show that the one-loop correction contains a hard 1/e UV divergence, as well 
as a logarithmic IR-divergent term of the type In \/iOp)^, thus considerably softening the 
usual UV/IR mixing phenomenon. We show that both of these problematic terms vanish 
for a certain choice of the noncommutative parameter 9 which preserves unitarity. We find 
non-perturbative modifications of the neutrino dispersion relations which are assymptoti- 
cally independent of the scale of noncommutativity in both the low and high energy limits 
and may allow superluminal propagation. Finally, we demonstrate how the prodigious 
freedom in Seiberg-Witten maps may be used to affect neutrino propagation in a profound 
way. 
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1. Introduction 

The study of spacetime quantization has originally been motivated by major problems 
of physics at extremely-high energies, in particular the problems of renormalization and 
quantum gravity. Heisenberg-type spacetime uncertainty relations can effectively lead to 
a replacement of the continuum of points by finite size spacetime cells, thus providing a 
means by which to tame UV divergences. A branch of mathematics arising from these 
motivations has come to be known as noncommutative geometry. It is reasonable to expect 
that noncommutative (NC) field theory models can provide some guidance for a deeper 
understanding of the structure of spacetime at extremely-high energies. In fact, these NC 
models appear quite naturally in string theory Q. The relevant scale of noncommutativity 
may very well be beyond direct experimental reach for the foreseeable future (except in 
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certain theories with large extra dimensions). Nevertheless, non-perturbative effects can 
nevertheless lead to profound observable consequences for low energy physics. A famous 
example is UV/IR mixing. Another striking example is the running of the coupling constant 
of noncommutative U(l) gauge theory in the simple star (-k) -product formalism [^. The beta 
function 

«»^) = -^f.^iv^ (1.1) 

of NC U(N) gauge theory is identical to that of ordinary SU(N) Yang-Mills theory for N > 
1, but in the noncommutative case it remains valid even for Abelian = 1 gauge theory. 
Hence the theory will suffer from asymptotic freedom 1^, |3|^. This result is manifestly 
independent of the scale of noncommutativity and thus remains valid even for vanishingly 
small (but non-zero) noncommutivity. 

In a simple model of NC spacetime local coordinates are promoted to hermitian 
operators satisfying spacetime noncommutative relations 

=ie^^, (1.2) 



where 6^^^ is real antisymmetric matrix of dimension length? . The commutator ( |1.2[ ) implies 
spacetime uncertainty relations 

Ax^Ax^ > ^ir^- (1-3) 

It is straight-forward to formulate field theories on such noncommutative spaces as a de- 
formation of the ordinary field theories. The noncommutative deformation is implemented 
by replacing the usual pointwise product of a pair of fields f{x) and g{x) by a ^-product 
in the action: 

f{x)g{x) {f*g){x) = f{x)g{x) + 0{e,df,dg). (1.4) 



The Moyal-Weyl ^-product is relevant for the case of a constant 0^*^ and is defined as 
follows: 

(/*5)(x) = /(x)et^^'"'^5(^). (1.5) 

(The ^-product has also an alternative integral formulation, making its non-local character 
more transparent.) The operator commutation relation ( |1.2| ) is then realized by the star(-k)- 
commutator 

[x^,x''] = [x^U''] = iO^''- (1.6) 

In analogy to the introduction of covariant derivatives in gauge field theory, the star- 
product can be promoted to a gauge-field dependent covariant star product. Together with 
a gauge-field dependent covariant integral measure this generally leads to a noncommuta- 
tive gauge field theory based on so-called Seiberg-Witten (SW) maps |l| . The resulting 
type of noncommutative quantum field theory has been studied for quite some time. 

In this construction the noncommutative fields are obtained via SW maps from the 
original commutative fields. It is important to note that there is typically some freedom 



^The negative /?- function is most significant in pure noncommutative gauge theory. When fermion fields 
are added the situation can change considerably Q. 
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in the choice of SW map and that there is no warranty that every change in the choice 
of SW map win lead to a physicahy equivalent theory: Deformation, like quantization, is 
usually not unique and different deformations can lead to physically inequivalent models. 
The deformed model is not uniquely fixed by its commutative classical (tree level) ^ — >■ 
limit. Different SW maps can behave like different quantization procedures. 

The perturbative quantization of noncommutative field theories was first proposed in 
a pioneering paper about fifteen years ago Q. Since then considerable efforts have been 
devoted to this subject. However, despite some significant progress like the models in 
m, 1^, ^, a complete understanding of quantum loop corrections still remains in general 
a challenging open question. This fact is particularly true for the models constructed 
by using Seiberg-Witten map expansion since the map was for a long time expressed as 
an approximative expansion in the noncommutative parameter ^, |lT], 12]. One 
loop quantum properties p|, M, pi, 111, P, 111, 0, M, |23, [23], as well as the 



related phenomenology |26, 27, 29, ^ 31 1 of the ^-expanded models, have also been 
investigated recently. 

The tree and one-loop analysis of the minimal NCSM truncated to first order in 
1 19, 20, 31 1, has shown that considerably different physical behavior can arise from differ- 
ent choices of the Seiberg-Witten map. Furthermore, an analysis of the photon two-point 
function in a ^-expanded model up to all orders via SW map |14]^ reveals that the pa- 
rameters that fix the choice of SW map are running coupling constants |15]. This again 
indicates physical differences among different SW map deformations. 

Some results about closed- form solutions and/or alternative ^-exact approaches, start- 
ing from exact solutions for the Seiberg-Witten map, have existed for quite some time |32 



33, 34, |35[]. Quite recently, ^-exact SW map expansions, in the framework of covariant non- 
commutative quantum gauge field theory |36|, were applied in loop computation [R^, IS^, 39 1 



and phenomenology |4C, 41 1. These more sophisticated theories differ quite drastically from 
their ^-expanded cousins, as they introduce in general a nonstandard denominator p9q into 
the loop integral. A few methods have been proposed so far to handle this problem: An 
expansion and re-summation with respect to 6 in the loop integral allowed some progress in 
1 37]; another approximative method which also looks promising is integration using para- 
metric derivatives. In general however, obtaining result in a closed form still remains a 



challenging problem ]38]. 

In this article we obtain a closed formula for the one-loop correction to the propaga- 
tor of a massless neutrino (neutral fermion) in the adjoint representation of U^(l) gauge 
group In the evaluation we combine parameterizations of Schwinger, Feynman, and 
a (modified) HQET parameterization ]^^, which was developed originally for heavy quark 
effective theory (HQET). This model is relatively easy to handle since both the gauge field 
and the fermion Seiberg-Witten map can be defined using generalized star products. The 
model is interesting in its own right for both theory and phenomenology: Its unexpanded 
version features a fermion-boson number symmetry which can cancel the leading order 



^To absorb the loop divergences, at each order in 9 the freedom in the choice of the Seiberg-Witten map 
has been used |l5| ]. 
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ultraviolet /infrared (UV/IR) mixing [^]. The model was also considered as an example 



for tree level neutrino-photon coupling via noncommutativity [11, 



The radiative corrections that we obtain contain in general both a hard UV term 
and a logarithmic IR singularity. At a special value of the noncommutative parameter 0, 
both singularities vanishes. Analyzing the poles of the resulting (finite) modified neutrino 
propagator reveals, depending on the energy regime, modes with either heavy masses or 
modes whose propagation depends on the preferred direction in space set by 6^'^. In the low- 
energy regime we find modes propagating superluminaly. These properties present some 
previously unknown features of 6-exact noncommutative quantum field theories (NCQFT). 

The article is structured as follows: In the following section we describe the actions 
of two alternative models which differ with regard to the choice of SW map and we give 
the relevant Feynman rules. Section 3 is devoted to the computation of the one-loop 
neutrino self-energy. Nonequivalent divergences and corresponding dispersion relations for 
both actions are described. Asymptotic dispersion relations are given for the low-energy 
and the high-energy regimes. Section 4 is devoted to discussion and conclusions. Relevant 
computational details of the nontrivial loop-integrals are given in two appendixes. 



2. Model 

In setting up our models we adhear to the following principles of 9-exact NCGFT: 

(i) The main principles that we are implementing in the construction of all of our 0-exact 
noncommutative models are: The standard field content and the commutative gauge sym- 
metry as the fundamental symmetry of the theory are fixed. 

(ii) In the construction of the noncommutative action, generically, electrically neutral mat- 
ter fields will be promoted via (hybrid) Seiberg-Witten maps to noncommutative fields 
that couple to photons and transform in the adjoint representation of U^(l). 

(iii) The different actions discussed in this paper are constructed by employing SW map 
freedom. 

(iv) The inclusion of all gauge covariant coupling terms into these actions is a prerequisite 
for reasonable UV behavior. 



2.1 Action 



Taking the above into account we arrive at the following model of the SW map type 
noncommutative U*(l) gauge theory with a gauge field coupled to a massless neutral 
fermion ^ via a star-commutator 12] 



(2.1) 



The action is defined in the usual way ]11|, |12|, |35 

1 



S 



(2.2) 
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with definitions of the noncommutative covariant derivative and field strength resembling 
the corresponding expressions of non-abelian Yang-Mills theory: 



D^^ = d^^ -i[A^^,^] and F^^ = d^^A^ - d^A^ - i[A^ * A^]. 

All the NC fields in this action are images of the corresponding commutative fields 
and ip under (hybrid) Seiberg-Witten maps. In the original work of Seiberg and Witten 
and majority of the subsequent applications, these maps are understood as power series 
of the noncommutativity parameter 9^'^. Physically, this corresponds to an expansion in 
momenta and is valid only for low energy phenomena. Here we give an alternative point of 
view and employ an expansion in formal powers of the gauge field and hence in powers 
of the coupling constant e. At each order in we shall determine 0-exact expressions. 
In the following we discuss the model construction for the massless case, and we shall set 
e = 1. To restore the coupling constant one simply substitutes by ea^ and then divides 
the gauge-field term in the Lagrangian by e^. 

Next step we expand the action in terms of the commutative fields and ip using the 
following SW map solution |Q 

1 

-( 
2 



A^= a^- ^O'Pa, *2 {dpa^ + fp^) + 0{a^) 



ap*2{dpa,y-k20atp) - [apdp%(j{duaa + fua) - dpOf^ipa^Ua-] > + 0{a^)ip. 



(2.3) 



Here the two generalized star products 



sm 



died2 



[f{x)gix)hix)] 



*3 



fix) -k2 g{x) = ""g^gg^ f{xi)g{x2) 

sin(M93)si^( 9ig02+a3) , 
2 2 



X\=X2=X 



+ {1 o 2} 



f{xi)g{x2)h{xz) 



(2.4) 
(2.5) 



are symmetric in their arguments, but nonassociative. The resulting expansion in the 
coupling constant defines the one-photon-two-fermion, two-photon-two-fermion and three- 
photon vertices ^-exactly. 

The expansion of the action is straightforward using the SW map ( |2.3|) . The first 
nontrivial contribution in the gauge field expansion leads to the following photon self- 
interaction terms 

Sg= j idpa, * K t + {0'"'ap *2 (d^a, + /,,)) * + 0{a^). (2.6) 
The photon-fermion interaction up to 2-photon-2-fermion terms is derived using the first 
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Figure 1: Three- and four-field vertices 



order gauge field and second order fermion field expansion 



-il^l^[a^'ie'^ai^2d,i>]--i>-i^'[e'^ai^2{dja^+fj^,)'li^]^^^^ 

+ ^-e'-^e^^ (^au *2 {dm + fu)) *2 dji^ + 2ai *2 {dj{ak *2 ditp)) - *2 (^fcO-,- *2 SiV^) 

+2ai-k2{dj{ak*2diip))-ai-k2{dkaj-k2di7p)+ [aidkip{djai + fji)-dkdiipajai]^^ 
+ V'C(a3)V'. 



(2.7) 



It is important to note that the actions (2.6) and (2.7) for gauge and matter fields respec- 
tively, are nonlocal objects due to the presence of the non-local (generalized) star prod- 
ucts -k, *2 and +3. The appearance of non- locality in the actions is expected to be reflected 
in corresponding quantum corrections to the neutrino propagator (2-point functions). 



2.2 Feynman rules 



From the interaction terms we read out the three vertices needed for one-loop two point- 
function computations |44| 
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sin ■ 



(2.8) 



+ ^72) [e>'''{q2q3q'i " '/I'/s'/sO + ^^^('Z2'Z3'Zr " qiq2q^3) + ^'^^(^I'Zs'Zs^ - q^q^) 

- g'^iqU + g?9t) - g^'iqlrs + ^i^T) - g^'iqlq^ + qM^) 



+ g3'(9t^3 + qiqi;) + q'i{qiq^2 + q^q'2) + ^l^t^Z^ + qX) 



(2.9) 



y3"'iqi,q2,p,p) = 2i{2 



^ ^P^7 - 2i 



—p 



sin 2^ sin 21^, , „ , sin 2i|£ sin 



+ / 



qi0q2 

p6p' .,■ qi6q2 



qi6pq29p' 



q20p' 

2+i>')p^p'' 



qi0q2p9p 
sinSlMsin^^ 
Pi9p'q29p 



j^^q^Ope^- -2ql^p''^ 



-2{q2-pW + 



sin sin 22fP 



^ sin2^sin2lM ^ Sin2l|22sin2|^ 



q20p'qi9p' 



q20p'p9p' 



qi9p 

{2p''q!^ + 9'"'p9q2 - p^p" 



+ i> 



+ 



sm ^-7^ sm 



i2q>^p''' - q29p'9n 



qi9q2p'9p 



^2iq2+pTP 



2 ^^^^ 2 _Qtiv 



qi9p 



qiopq2 

sin^sin23|£ sin gin £|£ \ , ^ ^ ^,„^ ,^ , , , ^ 
q29ppi9p q29ppt)p J 

+ {qi ^ q2 and jj. ■H' u} , 

where q^ = {9qY = 9'^''qu, and in addition we define q^ = {99qY = 9'^'^9^pqP . 

2.3 Alternative actions and Feynman rules 



(2.10) 



As discussed in [41], tliere exist alternative consistent and covariant choices for the non- 
commutative interactions. This is related to a freedom in the choice of SW maps used in 
the construction of the theory. For the action ( |2.7D and Feynman rules ( |2.8D , ( ^.10 ) , we 
presented in this section an alternative construction of the massless action, with coupling 
constant e = 1. 

We start with the action for a neutral massless free fermion field 



(2.11) 



Sf= / tp^'^d^tp = / V'^t'^s^V' 



where, as indicated, a Moyal-Weyl type star product can be inserted or removed by partial 
integration. Following the method of constructing a covariant NC gauge theory outlined 
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at the beginning of this section, we hft the factors in the action via (generahzed) Seiberg- 
Witten maps ^'[i/j, ...], ^[df^ip, ...] to noncommutative status as follows: 



Sf ^ = J ^WrHOf,^l^) = J ^W*j^'<^{d^,p) . (2.12) 

Now if the SW maps ^, $ and a corresponding map A for the gauge parameter A satisfy 

<5A(^W) = i[A(A)t^W], 
S^iHd.i^)) = ^[A(A) t , 

we will have a noncommutative action which is gauge invariant under infinitesimal commu- 
tative gauge transformations 5x and reduces to the free fermion action in the commutative 
hmit e ^0. 

The appropriate map ^ can be the same as in ( |2.3D . Recalling that we are dealing 
with neutral fields, i.e. dip = and 5{d^ip) = 0, we notice that we can in principle use the 
same map also for 

+ le^^-e^'|(afc *2 {dm + fu)) *2 djid^^P) + 2aMdj{ak*2di{d^m 

- ai*2(9fcai*29/(9^V')) - [akdi{df,7p){djai + fji) - dkdi{d^'4j)ajai] \ + Oia^)^) . 



(2.14) 



This construction is quite unusual from the point of gauge theory, as it yields a covariant 
derivative term without introducing a covariant derivative. In any case the resulting action 



y^e'^'liak M {dm + fu)) *2 dj^ + 2aMdj{akMdii^)) 



2 

ai*2{dkaj-k2ip) - [akdiip{djai + fji) - dkdi{'4>)ajai] (2.15) 



+ -0*^rV|(afc *2 {dm + hi)) *2 d,{^^) + 2aMd,{ak^2di{H))) 

- ai-k2{dkaj-k2di{0Tp)) - [akdi{0ip){djai + fji) - dkdi{^tp)ajai]^^ 
+ 0{a') , 

is consistent and gauge invariant. The action leads to the following photon- fermion inter- 
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action vertices, i.e. Feynman rule, 



Vr^Jq,p) = -F{q,p){eqr^, 
qi6q2pBp' 



(2.16) 



2f^sin2j£ _ sinai|l^sinaf^ \ 



— ^ 

g^ysin2f£ ^^, + p^^^ 

qiOp'q20p V 



sin ^^y- sm — 
qi0p'q26p 

/sin 2^ sin 2lM sin SiM gin 2^ \ ,^ „ 

sm 2 sin 2 



-f- 



sm^^V^sm^^V^ ^ /sm^^V^sm^ 

^ ^^2(^2 + p') V" - ' . ' 

V q20pqidp 



2 " 

qi9pq26p' 



■{2p^qi^ + rV^92 + P'^P I + {ft ^ 92 and/i 



There is also a second choice for 4>: 



(2.17) 



(2.18) 



This leads back to the original action discussed in the previous subsection. In general one 
can chose any superposition of the two SW maps ^ait2 and ^^am indicating a freedom in the 
choice of Seiberg-Witten map. Note that Feynman rule ( |2.8| ) is more natural from the point 
of view of gauge theory. In this article we analyze two different actions (Feynman rules), i.e. 
the original action and its alternative presented in this section. These essentially different 
actions ( |2.7D and (2.15), despite having the same field content and gauge symmetry, will 
be shown to produce different neutrino self-energies at one-loop in the next section. 

Finally, the pure gauge field action Sg in the alternative action, remains the same as in 
the original action. We keep Sg and ¥2'^'^ intact, that is we apply ( ^.91) in the computation 
of the two-point function from the alternative action. 



3. One-loop neutrino self-energy 
3.1 Diagrams 

According to the vertices defined in the last section, there are four diagrams which con- 
tribute to the neutrino self-energy at one-loop 



^l—loop 



Si + S2 + S3 + 



(3.1) 



The first Si and third S3 are the analog of the commutative fermion self-energy. The 
second one is similar to a tadpole graph in (j)'^ theory. However, it is not difficult to check, 
through direct computation from ( 2.10 ), that g^uV^'^ip, —p, q, q) = 0. Therefore S2 is equal 
to zero. The last two diagrams, S3 and S4 are tadpole diagrams obtained by connecting a 
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Figure 2: One-loop contributions to neutrino self-energy 



one-loop one-point function to a triple particle vertex. These two diagrams could be ruled 
out by the noncommutative charge conjugation symmetry defined in [^], so we do not take 
them into account. Namely, here we have taken the charge conjugation transformation to 

-e^"'. Thus only the 



be the same as in the equations (64) to (66) from |10|, i.e. 0'^^'^ 
two-point function bubble diagram Si needs to be evaluated. 



3.2 Two-point function 

The integral Si is defined as follows: 



4-D 



o '-^■F{q,p)[{qep)r+^q' 



(2vr) 

F{q,p) -{q9p)^P - i>f + if 



^{p + qf 



f d-q 




M 1 1 


1 {27t)D 


K 2^ ) 


g2 (p + Qr)2 1 



iq9p)\i-D)i^ + 



(3.2) 



+ (qOp) [mP^ + 2p-q)- 0{2q^ + 2p ■ q)] 

+ [Afip"^ + 2p-q)- q^{f + 2p ■ q)) + i{f{q^ + 2p ■ q) - p^{f + 2p ■ q) 



where n represents a dimensionful regularization parameter. Note that during the evalu- 
ation of Si we shall hide the explicit e-dependence in the equations. We shall later make 
the e-dependence explicit during the derivation of the dispersion relations. 

Combining the power of the (qOp) term and the phase factor from F"^, we can separate 
Si into a sum of three integrals 



-Si = {4-D) /i23 + /456 + /■ 



789 



(3.3) 
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-^123 
-^456 
-^789 



iqOp ^—iqOp 



1 



1 



d^q 1 1 



(3.4) 

[|(2p2 + 2p-q)- pi2q^ + 2p ■ q)](2 - e^^'P - e-*'?^^'),(3.5) 



{p' + 2p-q)-q\p^ + 2p-q)) 



{2tt)D q2(p + q)2 ^q0p)2 
+ ^ip\q'^ + 2p-q)- p\f + 2p ■ q)) 



2-e 



(3.6) 



where the presence of the regularization mass parameter fi in the above integrals is un- 
derstood. The exphcit evaluation of the above loop-integrals, as given in Appendix A, 
yields 



Si 



Ytr6'6' 2p 



, o / D\ [D D 



2 j dx{l - x) {x{l - x)/) 4 ^ 2^-'^{p^f-^K^_D (x(l - x)p^f) 



(4- 



TT 



2sin^ 



2\2- 



dx{\ — x){p ) 



(x(l — x)p p ) 



IF2 



2 1" ^ V 2' 2' 2 ' 



1 3 L> x(l — x)p^fP' 



iD-2t 



3-Z? 



1-^2 



3-Z) 4-1) 5-1) x(l- x)pV 



(3.7) 



3.3 Neutrino self-energy 



The general expression for the neutrino self-energy, after evaluating integrals in (| 
D = A — e and in e ^ limit, receives the following closed form structure 



for 



Si = -7m 



p^ A + {eopY 



(Op)' 



B 



(3.8) 
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with 



A 



1 



(47r)2 
1 

(4^ 



(Sl + 2s2)Ai + (1 + Si + S2)A2 
>1 + 2S2 



(3.9) 



^ln^ + 2-,oQ)+ln(4.) 



^ 2^rfA: + l 



1 

7r2 



-r(i,- + i)r(fc + |) 



+ 



(4vr) 



KI + S1 + S2) 



16 ; V 16 

2^^r(fc + i)r(fc + f)r(fc + |) 



16 



2\ 



In 



V'o (fc + 1) + V'o ( A; + M - f A; + ^ j - Vo f A; + ^ 



16" 
(3.10) 



-1 



(4vr) 



; V vrir(fc + i) 



16 



2\ 



In 



V'o (fc + 1) + -(Ao U + 



1 



^0 u + 



2J ' V 2 
Here si^2 are scale-independent ^-ratios 

2 tr^e 

si=p 



16 



r^2 



(4^) 



S2=P 



(3.11) 



(3.12) 



(0p)2' ^ (0p)4 

and 2 correspond to the first and second square brackets in Eq. ( 3.10 ), respectively. 

It is important to note here that amongst other terms contained in both coefficients 
Ai and A2, there are structures proportional to 

(/(M^)"^^(ln(p^(M^))"', Vnandm = 0,l. (3.13) 

The numerical factors in front of the above structures are rapidly-decaying, thus the series 
are always convergent for finite arguments as we numerically demonstrated below: 



Ai ~ - + \n{^Te^^^?{epf) 
11 



(3.14) 



7129 



p2(0 )2 f ^ 2(^ )2) ^^( 2(^ )2)2 

72^ ^ \ 8800^ ^ ' ' 313600^ ^ ' ' 17882726400 



{p\epYY + ... 



12 



80 



13440 



3870720 



A2 = -Svr^S 

^ 2 + -p^epf ( 1 + -^pHopf + {pHopff + — ^ — 

18 ^ ^ V 8400 ^ ^ 21952000^ ^ ' ' 19914854400 



(3.15) 



{p\Bpff + 



1 1 + In ^\ p'^^p^' (i + p'(^p)' + (p'(M')' ^ ipH0p?r ^ 



4 y / 12 V 22400 6773760 

where 7e — 0.577216 is Euler's constant. It is to be noted here that the spinor structure 



proportional to /J is missing in the final result, confirming thus conclusion from |45|. 
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3.3.1 Divergences and counter terms 

The first striking fact of our closed form result is the existence of a non-local UV divergence 



term 



^UV 



p 



tree ^^{eepf 



(4vr) 



(3.16) 



XepY {epY 

which does not vanish in the 0—7-0 limit. This term also clearly differs with respect to 
a model not based on Seiberg-Witten maps, where the UV divergence does not have the 
momentum and e depended factor (si + 2s2). Therefore, the existence of such a divergence 
suggests a necessity to introduce the following nonlocal counter-term 

{eedf' 



162 



tiee 



+ 2- 



{edy (ed) 



(3.17) 



which would cancel it. 

Besides the hard 1/e UV divergence, there is a soft UV/IR mixing term 

{eepf' 



^UV/IR = 



p 



tree 



(4vr) 



In 



16 



(3.18) 



(6*^)2 (6*^)4 

represented by a logarithm, and it diverges at both the ultraviolet and infrared limits. 
Since the soft UV/IR mixing ( p. 181 ) appears in ( p. 10 ) with exactly the same coefficient as 
the UV 1/e term does, the introduction of the same nonlocal counter-term, that is ( |3.17] ), 
would remove it. 

Finally, both these terms, ( ^.16 ) and ( p.l8| ) respectively, are proportional top^. There- 
fore if the counter term ( p. 17 ) is included and the renormalization point is selected at |^ = 0, 
our result indicates that the dispersion relation p^ = could still hold. However, in the 
next subsection we shall investigate the other solutions too. 

In the renormalization procedure, all three coefficients Ai, A2 and B from ( p.lOj ) and 
(3.11) are renormalized by subtracting counter terms from Si in (p.^). We obtain then the 
renormalized neutrino self energy 

p^ Ar + {eepr 7^ Br 



'IR 



Ar 



1 



(47r)^ 



(0p)2 J 
(.Si + 2S2)AiR + (1 + Sl + S2)A2R 



(3.19) 
(3.20) 



1,2; 



(3.21) 
(3.22) 



where 

ARip^ {ep)\ {eepf) = A{p^ {ep)\ (eepf) - a,{pI {ep)l {eep)l) 
Br{p\ {ep)\ {eepf) = B{p\ {epf, {eepf) - b{pI {ep)l {eep)l) , 

and (pq, (6*^)0, {eep)^) is a choice of renormalization point. 

For unitarity [^] reasons, it is convenient to take the following choice of degenerate e 

/O o\ 
10 
0-100 
\0 0/ 



1 



(3.23) 
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One finds that si + 2s2 = in this degenerate case, therefore Suv = 5] 



counter term is needed. However, note that in this case (Op)"^ 
renders the modified propagator zero in this subspace. 



UV/IR 
for pi = p2 



■ 0, so no 
0, which 



3.3.2 Dispersion relation 

In order to probe possible physical consequence of the one-loop quantum correction S]^_joop, 
we consider the modified propagator 

1 



1 



l-loopM 



S2 ■ 



(3.24) 



The Minkowski counterpart Si-Zoop^^ of (|3.lD is computed by using the established param- 
eterization technique and Wick rotation. The resulting one loop correction is 

n2 



^1-loopM 



p^" A + {eopf 



p 



{OpY 



B 



(3.25) 



We further choose the noncommutative parameter to be ( p. 231) so that the denominator of 



(3.24) is finite and can be expressed explicitly: 

E2 



P 



^2 



4 + 24 + 5 

Pr Pr 



^2 6 + 2 



+ 1 



(3.26) 



where pr represents r-component of the momentum p in a cylindrical spatial coordinate 
system and A2 = e'^ A2 / {4:Tr)'^ = -B/2. 

We are interested in the zeros of the denominator, especially the simple zeros which 
have the local form pQ — f{pi,9), for this is associated with the time-evolution of the 
corresponding excitation via the Fourier transformation of the propagator 



Y[ dpi J dpo ^ e 



ipo{x—x')^—ipi(x—x'y 



(3.27) 



In general, a factor exp[if{pi,9)t] will arise from through the residue at the zero-point 
Po = f which stays in the upper half of the complex plane. The existence of the nonzero 
lm[/] is also natural since it implies that this excitation decays with respect to time and 
thus, it is unstable (the tachyonic modes with (pq + m^) type pole can be considered as a 
special case). 

From above one see that p^ = defines one set of the dispersion relation, corresponding 
to the dispersion for the massless neutrino mode, however the denominator ( p. 26 ) has now 
one more coefficient 

^2^ 



^2 



4 2 
^ + 2^ + 5 

Pr Pr 



A2 6 + 2 



+ 1, 



(3.28) 



which could also induce certain zero-points. Since the A2 is a function of a single variable 
p'^Pr, with p2 = Po — Pi — P2 ~ P3 and p^ = Pi + pIj the condition S' = which we are 
interested in can be expressed as a simple algebraic equation 



Alz' 



2 Ao 



^2 



z+ (l-6A2 + 5Ai 



0, 



(3.29) 
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of new variables z '■= /pi, in which the coefficients are all functions oi y : = p^p^/hl^ 



NC- 



1 - As ± 2 



(3.30) 



The two formal solutions of the equation (3.29) 

1 

are direction dependent, i.e. birefringent. The behavior of the birefringent solutions ( 3.3C| ), 
with respect to a propagating energy, can be analyzed at two limits y — > 0, and y — >■ oo. 



The low-energy regime: p^p^ <^ At 



NC 



For y ^ 1 we simply set A2 and A2 to their zeroth order value e /Svr and e /647r , then 



8^^ 



1 ±2 



8^^ 



(3.31) 



or equivalently 




1 ±2 



8^^ 



1 • p,? ~ (859 ± 59) • p"}. , 



(3.32) 



defines two (approximate) zero points. Prom the definition of p"^ and we see that both 
solutions are real and positive. Taking into account the higher order (in y) correction these 
poles will locate nearby the real axis of the complex po plane thus correspond to some 
metastable modes with the above defined dispersion relations. As we can see, the modified 
dispersion relation ( |3.32| ) does not depend on the noncommutative scale, therefore it intro- 
duces a discontinuity in the Anc ~^ 00 limit, which is not unfamiliar in noncommutative 
theories. 



The high-energy regime: p'^p'^ » A^q 

At y ^ 1 we analyze the asymptotic behavior of A2 starting with its integral form 



TT 



2 sin |7r 



A'' 



dx{l 



x{l — x)p p. 



.2^2 \ 1 



22-ep 



Pr 



4 + 



x(l — x)p'^p^ 

.2^2 



1 eel e x{l — x)p p: 
'2 ^2' 2' 2 ^2' 4A^ 



"-NC 



A^ 



3 

■7T2 



p^pI 

16A^c 



1- 



2 fl 

r . 

2 



1F2 



13 5 



1 e 

■2 + 2 



2-^3 



"2'^2' '2'2"^2'2' 



p^pI 



p^pI 

2' 2' 2 ~ 2'~16A^ 



Using the documented asymptotic expansion of generalized hyper-geometric series 
leading and next-to-leading asymptotic orders of A2 when y — >■ 00 reads 



(3.33) 
1, the 



A, 



i-K^ 1 



1 — IQin y e 



+ 0(y 



(3.34) 
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So at leading asymptotic order 

z 1 ± 2i, 



(3.35) 



or 

vl^vi±2ipl. (3.36) 

We thus reach two unstable deformed modes besides the usual mode = in the high 
energy regime. Here again the leading order deformed dispersion relation does not depend 
on the noncommutative scale Anc- 

3.4 Neutrino self-energy for alternative action 

Using the Feynman rule ( |2.16| ) of the alternative action, we find the following closed form 
contribution to the neutrino self-energy from diagram Si 



Si 



(4vr) 



+ 4- 



3(ep)2V(M^ {Qvf 



(3.37) 



The detailed computation is presented in Appendix B. From ( p.l7| ) one gets the relation, 
9^^v^z^^^J<P■,~'P^(l^^) = 0> showing that S2 = 0, while S3 and S4 vanish due to charge 
conjugation symmetry. Therefore we have again Si_/oopa;t2 ~ ^iaH2- There is no alternative 
dispersion relation in degenerate case ( ^.23 ), since the factor that multiplies |5 in ( 3.37] ), 



does not dependent on the time- like component po (energy). 

We have to notice that equation ( |3.37| ) is much simpler than the corresponding expres- 



sion (3.8). This is not unfamiliar for actions arising from different SW map deformations. 
There are no hard 1/e UV divergent and no logarithmic UV/IR mixing terms, and the 
finite terms like in A\ and A2 are also absent. Thus the subgraph Si does not require any 
counter-term. However, the result of the subgraph Si evaluation, from alternative action 2 
(2.15), does express powerful UV/IR mixing effect due to scale dependent ^-ratios. Namely, 



in terms of scales only, the Si^^^^ experience the forth-power of the NC-scale/momentum- 
scale ratios ~ in ( 3.37| ), i.e. we are dealing with the Si^^^j ^ ^i-^Nc/p) within 



the ultraviolet and infrared limits for Anc and p, respectively. 

The absence of new spinor structure in the alternative neutrino self-energy (3.37) fur- 
ther suggests the possibility of an appropriate field strength renormalization with suitable 
divergence cancellation for 9^0 limit. Here certain hint may be found in the counter 
terms proposed for the translation invariant renormalizable noncommutative (f)'^ model with 
regular commutative limit H. 



4. Discussion and conclusion 

In this article we present a 0-exact evaluation of the one-loop quantum correction to the 
neutral fermion propagator. We in particular evaluate the neutrino two-point function in 
a two different Seiberg-Witten map based NCj^QED models. As there is no a priory way 
to exclude either model, we study them both and point out their different behavior. Our 
method, based on 0-exact expressions of Seiberg-Witten maps, together with a combina- 
tion of Schwinger, Feynman, and HQET parameterization, yields the one-loop quantum 
correction in a closed form for both models. 
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From the neutrino one-loop two-point function we obtain the self-energy of a massless 
neutrino and a dispersion relations which depends on spacetime noncommutativity. The 
general expression for the neutrino self-energy given in ( |3.8[ ) contains both a hard 1/e ul- 
traviolet term ( |3.16D and the celebrated UV/IR mixing term with a logarithmic infrared 
singularity In \J (Op)"^. The later reflects the fact that the UV divergence is at most logarith- 
mic, i.e. there is a soft ultraviolet /infrared term representing UV/IR mixing ( 3.18| ), which 
is similar to the logarithmic term in the usual vacuum polarization of the photon, in simple 



^-product based NCQED without SW maps [48|. Since we have already discussed in detail 



properties of the UV/IR mixing term ( 3.18 ) in [39| we shall not repeat that discussion here. 



The essential difference of our results (3^) as compared to |48, 49, 51, 52 1 is that 



in our case both terms are proportional to the spacetime noncommutativity dependent 
^-ratio factor (si + 2^2), which arise from the natural non-locality of our actions and does 
not depend on the noncommutative scale, but only on the scale- independent structure 
of the noncommutative 0-ratios. This behavior, being non-perturbative in nature, differs 
from that of fermions in the fundamental or in the adjoint representation in usual, 
product only based and 0-unexpanded NCQED. Besides the divergent terms, a new spinor 
structure {66p) with finite coefficients emerges in our computation, see ( |3.8| )-( |3rT^ ). All 
these structures are proportional to p^, therefore if appropriate renormalization conditions 
are imposed, the commutative dispersion relation = can still hold, as a part of the full 
set of solutions obtained in (|3.26| ). Some of the propagating neutrino modes acquire mass. 
As the mass depends on the direction of propagation with respect to the noncommutative 
background set by 9^'^, these modes are birefringent, thus confirming previous result for 
chiral fermions in NCQED at first order in 6 ||25[| . 

The alternative action ( p. 151) , presented in section 2.3, has the same field content 
and gauge symmetry as the action ( |2.7| ), but a different choice of deformation freedom. 
Consequently these two different 0-exact actions led to two different neutrino self-energies 
( |3.8| ) and ( 3.37 ), respectively. For the unitary choice ( |3.23| ), and in the limit — > 0, 
self-energy (^]^) is finite, however the alternative one, ( 3.37] ), clearly diverge. This fact 
indicates strongly that the above two actions are not related by a field redefinition^. The 
corresponding alternative neutrino self-energy (|3.37| ), has less striking features than (^), 
but it does have it's own advantages owing to the absence of a hard UV divergences and the 
absence of complicated finite terms. Also, there is no modification of neutrino dispersion 
relation from ([3.37] ) in degenerate case ( |3.23| ). The structure in ( 3.37| ) is different (it is 
NC-scale/ energy dependent) with respect to the scale- independent structure (si + 2^2) 
from ( |3.9| ), as well as to the structure arising from fermion self-energy computation in the 
case of ★-product only unexpanded theories [|^. However, ( ^.37 ) does posses powerful 
UV/IR mixing effect. This is fortunate with regard to the use of low-energy NCQFT as 
an important window to quantum gravity |53] and holography ]5^ . 

The low energy dispersion relation (3.32) is, in principle, capable of generating a direc- 



'^It was shown for NCQED that only at first order in 9 the Seiberg-Witten map is a field redefinition, 
while at higher order in 9 the SW map can not be regarded as the field redefinition jIHj. This is certainly 
also true for the SW 6'-exact models. For S-exact models different SW maps are very much like different 
quantization procedures. See discussion of this issue in the Introduction. 
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tion dependent superluminal velocity, this can be seen clearly from the maximal attainable 
velocity of the neutrinos 



where i? is the angle with respect to the direction perpendicular to the NC plane. This 
gives one more example how such spontaneous breaking of Lorentz symmetry (via the 6- 
background) could affect the particle kinematics through quantum corrections (even with- 
out divergent behavior like UV/IR mixing). On the other hand one can also see that the 
magnitude of superluminosity is in general very large in our model, thus seems contradict- 
ing various observations. ^ The authors consider here, on the other hand, that the large 
superluminal velocity issue may be reduced/removed by taking into account several further 
considerations: 

• The reason to select a constant nonzero 9 background in this paper is computational 
simplicity. The results will, however, still hold for a NC background that is vary- 
ing sufficiently slowly with respect to the scale of noncommutativity. There is no 
physics reason to expect ^ to be a globally constant background ether. In fact, if 
the 9 background is only nonzero in tiny regions (NC bubbles) the effects of the 
modified dispersion relation will be suppressed macroscopically. Certainly a better 
understanding of possible sources of NC is needed. 

• In our computation we considered only the purely noncommutative neutrino-photon 
coupling, it has been pointed out that modified neutrino dispersion relation could 
open decay channels within the commutative standard model framework [^^. In 
our case this would further provide decay channel(s) which can bring superluminal 
neutrinos to normal ones. 

• Finally, as we have stated in the introduction, model 1 is not the only allowed 
deformed model with noncommutative neutrino-photon coupling. And as we have 
shown for our model 2, there could be no modified dispersion relation(s) for defor- 
mation(s) other than 1, therefore it is reasonable to conjecture that Seiberg-Witten 
map freedom may also serve as one possible remedy to this issue. 
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A. Loop integrals 



A.l Integral I123 



Integral I123 follow the same computation in NCQED without SW map |48], resulting in 



'123 



-2|i (p2)T-2r 2-- 



(47r)2 V 
1 

2 j dx{l-x) [x{l-x)p^)' 




2 



(x(l — x)p p ) 



(A.l) 



A. 2 Integral parameterizations for 1456 and lygg 

Next part involves integrals 1456 and Jygg which differ from I123 by the existence of a 
non-quadratic qOp denominators. To overcome this problem we introduce the HQET pa- 
rameterization |B3], represented as follows 



1 



ni ?12 



r(ni)r(n2) io (iaiy + a2)"^+"2 



(A.2) 



To perform computations of integrals (3.4), ( |3.5| ) and ( |3.6D , we first use the Feynman 
parameterization on the quadratic denominators, then the HQET parameterization help 
us to combine the quadratic and linear denominators. So for the denominators in I456 we 
have 



1 00 

II./'/" 1 
^ 2^ / dec / dy 

q^ip + qfqOp J J [(q2 + i^^f^i _ x) + {{p + q)'^ + ie)x + iy{qep)f 



2i dx dy 



1 



{q + xp + iyp)"^ + x{l — x)p^ + ^p^ + ie 



(A.3) 



Now we use the Schwinger parameterization to turn the denominators into Gaussian 
integrals. This then combines with different phase factors. For the zero phase, after an 
integral over the loop-momenta I = q + xp + ^yp and changing variable y to y = ay we 
have 



1 00 
I4 = 2i J dx J dy 



1 00 



72 



• exp 



d^l 



P{%-2)p-iy{l-x)pp^ + ^) 



^2 



D 



(2^) 



d^l 

D 



P + x{l — x)p^ + + ie 



(2vr) 



daa 







y2 

a/^ — ax(l — x)p'^ — : — p^ + ie 
4a 



(A.4) 
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For the other two nonzero phase factors, one can first performe the Schwinger parameteri- 
zation Hke in the zero-phase case, then change variable y to F = a (y =b The resulting 
integrals are 



'5.5±0.5 



1 oo 

= i I dx I dY 



±1 



(2vr, 



daa 



2 



2a2^ 



exp 



y2 

—ai^ — ax{\. — — — i?^ + 

4a 



(A.5) 



One can easily see that they differ from the zero-phase case just by a different Y integral 
limit. Summing over these three terms and integrating over the loop-momenta, one obtains 



-^456 = 2/4 — (is -I- /e) 

1 1 

= 2 — I dx dY / daa~^Y{l — x)pp"^ ■ exp 



(47r)f 



1 ipp^ 
(47r)f P'^ 



1 00 



000 

1 00 



y2 ■ 

—ax(l — x)p^ — - — 
4a 



{/dx/daa-#(l-x)exp[-ax(l-xy 







-/dx/A..'-f(l-x)exp 







— ax(l — x)p^ 

4a 



(A.6) 



The integral over a is then straightforward, so we have 



-^456 



1 4^p^ 
(47r)f P'^ 



1 00 



I dx daa^ 2 (1 — a;) exp [— ax(l — x)p^] 







1 00 
d./daa-f(l-x)exp 


1 4pp'^ ( 2 



-axil — xjp — — 
4a 



D D 



(A.7) 



(4.)f P- i^^^^^-^H^-2jn2-s 

1 

2 J dx{l-x){x{l-x)p'^)^~^ ■2^-^{fy-TK^_n |. 



A. 3 Integral I^sg 



A parameterization procedure introduced in the previous subsection yields the following 
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outcome for the integral Ijsg- 



'789 



D 



(2vr) 



dx / dY / daa exp 





2 



-a 



+ ie) + x{l - x)/ + 



;2m 



P(l - x)y ( ( ^ - 1 j - 



-y3(i-x)|(pV-i^') 



(A.8) 



The rest of the integrals are much more complicated but can still be handled. The integral 
/ygg after the Gaussian integration over loop- momenta I reads as follows: 



-^789 — — 



(47r)# 



dx / dY / da exp 



-a ( x{l-x)p^ + Y^'^ 



5^2 M 

4 



(A.9) 











(1 - X)ya2-f ( 1 _ + p^f - A - (1 - X)y3«3-f ^ (^2^2 _ 



+ p2| _ ^p2^^ _ (1 _ x)y2«2-f ^ (^V -P') } • 



To evaluate ( A.9 ), one first has to separate it into two parts 

1 1 + T 



(A.IO) 



where 



-2 j dx j dY j daexp -a (^x{l - x)p'^ + Y'^^^ (A.ll) 



(1 - x)ya2-f (^2^ P'P - ^-^P^A - (1 - x)y3a3-f ^ (^V - p") 



I' = Xfya^-f ^ ai-f'l . 



(A.12) 



The integral X can be evaluated in the same way as 1456 since 



I 



1 

1 ^ 00 
2 





(l-x)a2-f UJl-^pj +/|_^^p2 A _(^_^)^2^3-f^(^2^2_ 



(A.13) 
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The resulting Y and a integration are given in terms of the integrals over the modified 
Bessel functions K^^dIx] and K-^^_d[x], respectively: 



I 



Ml 



2 f- 2 p2 



J - ^ (P P -f>) 



(A.14) 



2 ' 2 



-2 / dx{\ - x) {x{\ - x)p^) 4 ^ 2^~'^{f)^~^K^_ii \[x(\ - ' 



-2^ (^V _ ^ ^^^^ _ ^^^^ _ ^^)p2^ ^ 



2f-i(p2)|-f [(x(I-x)pV) 



Next we evaluate the integral T' from ( [A.12[ ), where one could again first integrate over 
y, which yields both, the error function erf[x] and the modified Bessel function i^i 

2 

respectively 



tvoe 



(A.15) 



■{p'^) 2 / dx{l — x) / daa^ 2 ■ exp [— a — erf 



2ah 



+2^(^V-P^) / dx{\ - x) {x{\ - x)p^Y ^ 2^~^[f)^-^K-^_D [(x(I-x)pV) 
p J 2 L 



The integral over the error function can be expressed in terms of a regularized hyper- 
geometric functions 



vra / dx{l — x) j daa'^ 2 • exp [— a (x(l — • erf 



2a5 



2 sin 



dx{l — x)(p^)2 • 



--1 



^'2 



1^^ [I /i D x(l -x)pV 
2' 2'T' 4 



D-2/--2\l-4f 



-r 



(A.16) 
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Therefore, in the end we have 
1 



-^789 



D 

(47r) 2 



, y,D ^ I D\ f D D 



- 2 / - x) (x(l - x^) 4 ^ 2T-2(p2)i-fi^2-D - x)pV) 



)-^(*V-p')} 



vr 



2sin^ 



i 

j dx(l-x)(p2)2-f 



(x(i - ' r - 1F2' • • ^ ^ 



_2^-2r ('^hi^2 



2y ^ " V2' 2' 2 ' 4 
3-D 4-D 5-i5 x{l-x)p^p^ 



(A.17) 



B. Two-point function of alternative model 

In this appendix we compute the one-loop integral Si of the alternative model 2, which 
starts with 



D„ x2 sm 



2 / qOp 



By the parameterization we split S into the two integrals, I and I' , respectively: 



(4vr: 



= 1-1', 



(B.l) 



(B.2) 



^2 ^ y2|2 

4a2 





The Gaussian integration over / than produces 



1 



(47r: 



J dY-Y J da 

' ^ 



/tide 



a 2 -I ^ 2 ) e 4a 



/' = I{dY -Y ^dY) 



(B.3) 



(B.4) 
(B.5) 
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The change of coordinate a 
2 



(4vr: 



tvee f 

4 



1/a and then the integral over p gives 

-, 1 



D 



p2 
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/' = i{dY • y^--^ ^ dy • y-^) 



Now we assume 



dY ■ yi-^ 



D 



dY ■ Y-^ 



1 - D 



(B.6) 
(B.7) 

(B.8) 



This is clearly incorrect for D = 4 as the integral diverges at the lower limit. Here our 
motivation is from the general practice of dimensional regularization, in which the integrals 
could be evaluated at a certain D value where it converges, then take the extended D value. 
We then set D = A and finally obtain 



-1 



(4^) 



I' 



-I. 



(B.9) 
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